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Abstract. All existing "positive" results on two neutrino double beta decay in different nuclei were 
analyzed. Using the procedure recommended by the Particle Data Group, weighted average values 
for half-lives of '**^Ca, ^^Ge, ^^^Se, '^'^Zr, ^^^Mo, ''^"Mo - ^^^Ru (Of), "^Cd, '^"Te, ^^ONd, ""^^m - 
sm (0+) and 238u were obtained. Existing geochemical data were analyzed and recommended 
values for half-lives of '^^Te, '^"Te and '^"Ba are proposed. We recommend the use of these results 
as presently the most precise and reliable values forhalf-lives. 
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INTRODUCTION 

At present, the two neutrino double beta (2v/3/3) decay process has been detected in a 
total of 10 different nuclei. In ^*^^Mo and ^^°Nd, this type of decay was also detected 
for the transition to the 0+ excited state of the daughter nucleus. For the case of 
the ^^*^Ba nucleus, evidence for the two neutrino double electron capture process was 
observed via a geochemical experiment. All of these results were obtained in a few 
tens of geochemical experiments, more then thirty direct (counting) experiments, and 
in one radiochemical experiment. In direct experiments, for some nuclei there are as 
many as seven independent positive results (e.g., ^*^^Mo). In some experiments, the 
statistical error does not always play the primary role in overall half-life uncertainties. 
For example, the NEMO-3 experiment with ^^*^Mo detected more than 219000 useful 
events [T], which results in a value for the statistical error of ~ 0.2% . At the same time, 
the systematic error in many other experiments on 2v/3/3 decay generally remains quite 
high (~ 10 — 30%) and very often cannot be determined very reliably. As a result, it 
is frequently quite difficult for the "user" to select the "best" half-life value among all 
existing results. In fact, however, using an averaging procedure, one can produce reliable 
and accurate half-life values for each isotope. 

In the present work, a critical analysis of all "positive" experimental results has been 
performed, and averaged (or recommended) values for all isotopes have been obtained. 

The first time that this type of work was done was in 2001, and the results were pre- 
sented at MEDEX'Ol [2]. Then upgrated half-life values were presented at MEDEX'05 
[111]. In the present paper, new positive results obtained since 2005 have been added and 
analyzed. 



PRESENT EXPERIMENTAL DATA 



Experimental results on 2v/3/3 decay in different nuclei are presented in Table 1. For 
direct experiments, the number of useful events and the signal-to-background ratio are 
presented. 



DATA ANALYSIS 

To obtain an average of the available data, a standard weighted least-squares procedure. 



as recommended by the Particle Data Group [1450, was used. The weighted average and 



the corresponding error were calculated, as follows: 

X±dx = £w,.^,/£w, ± (£w,)-l/2^ (1) 

where Wi = l/ (dxi)^. Here, Xj and dxi are, respectively, the value and error reported by 
the i-th experiment, and the summations run over the N experiments. 

The following step is to calculate = ^^d compare it with N - 1, which 

is the expectation value of if the measurements are from a Gaussian distribution. If 

/{N —I) is less than or equal to 1, and there are no known problems with the data, 
we accept the results. If x^ / {N —I) is very large, we may choose not to use the average 
at all. Alternatively, we may quote the calculated average, while making an educated 
guess of the error, using a conservative estimate designed to take into account known 
problems with the data. Finally, if X^/{^ ~ 1) is larger than 1 but not greatly so, we 
may still average the data, but can increase the quoted error, 5x in Equation 1, by a scale 
factor S defined as 

S=[xViN-l)V^^. (2) 

For averages, we add the statistical and systematic errors in quadrature and use this 
combined error as 5xi. In some cases only the results obtained with high enough signal- 
to-background ratio were used. 



There are three independent experiments in which 2v/3/3 decay of ^^Ca was observed 
The results are in good agreement. The weighted average value is: 

ri/2 = 4.4^gJ.ioiV 



76Ge 

Let us consider the results of five experiments. First of all, however, a few additional 
comments are necessary: 

1) Result of the Heidelberg -Moscow group was corrected. Instead of the previously 
published value T1/2 = [1.55 ±0.01{stat)tQ \lisyst)] ■ lO^i y [|46i], a new value T1/2 = 



[lJ4±0.0l{stat)tQ\lisyst)] ■ 10^^ y lljj] has been presented. It is the latter value that 



has been used in our present analysis. At the same time, using an independent analysis. 



the Moscow part of the Collaboration obtained a value similar to the result of Ref. [|12ll . 
namely T1/2 = [lJS±O.Ol{stat)tQfoisyst)] ■ lO^i y 

2) In Ref. [H, the value T1/2 = O.Ql+J JJ • 10 y was presented. However, after a more 

careful analysis, this result has been changed to a value of T1/2 = I-^Iqi ■ 10^^ y llioll . 
which was used in our analysis. 

3) The results presented in Ref. [7] do not agree with the more recent and more 
precise experiments Furthermore, the error presented in [0] appears to be too 



small, especially taking into account the fact that the signal-to-background ratio in this 
experiment is equal to ~ 1/ 10. It has been mentioned before [48] that the half-life value 
in this work can be ~ 1.5 — 2 times higher because the thickness of the dead layer in 
the Ge(Li) detectors used can be different for crystals made from enriched Ge, rather 
than natural Ge. With no uniformity of the external background, this effect can have an 
appreciable influence on the final result. 

Finally, in calculating the average, only the results of experiments with signal-to- 



background ratios greater than 1 were used (i.e., the results of Refs. [|12l. lid. Illll). The 
weighted average value is: 



ri/2 = (i.5±o.i)-io2V 



82Se 

There are three independent counting experiments and many geochemical measure- 
ments (~ 20). The geochemical data are neither in good agreement with each other nor 
in good agreement with the data from direct measurements. Formally, the accuracy of 
geochemical measurements is typically on the level of a few percent and sometimes 
even better. Nevertheless, the possibility of existing large systematic errors cannot be 
excluded (see discussion in Ref. [49]). It is mentioned in Ref. fsO] that if the weak inter- 
action constant Gf is time-dependent, then the half-life values obtained in geochemical 
experiments will depend on the age of the samples. Thus, to obtain a "p resent" half-life 
value for ^^Se, only the results of the direct measurements l4\ were used. The 



result of Ref. 115 IH is the preliminary result of [14], hence it has not been used in our 
analysis. It is interesting to note that the "lower" error in Ref. [14] appears to be too 
small. Indeed, it is even smaller than the statistical error, and that is why we use here 
a more realistic value of 15% as an estimation of this error. As a result, the weighted 
average value is: 

ri/2 = (0.92 ± 0.07) -lO^V 



96 



Zr 



There are two "positive" geochemical results [[itI . [18|| and two results from direct 
NEMO-2 [16] and NEMO-3 [Q] experiments. Taking into account the comment in 



section 3.3, we use the values from Refs. lfl6l [6ll to obtain a "present" weighted half- 
life value for ^^Zr of: 

ri/2 = (2.3±0.2)-10iV 



ioomo 

Formally, there are seven positive result^ from direct experiments and one recent 
result from a geochemical experiment. However, we do not consider the preliminary 
result of M. Moe et al. [20] and instead use their final result f2^, plus we do not 
use the geochemical result (again, see comment in section 3.3). Finally, in calculating 
the average, only the results of experiments with signal-to-background ratios greater 



than 1 were used (i.e., the results of Refs. 
used the corrected half-life value from Ref. 
decreased by 15% because the calculated efficiency (by MC) was overestimated (see 
Ref. [1530). Secondly, the half-life value was decreased by 10% taking into account that, 
for the special case of ^^*^Mo we have to deal with the Single State Dominance (SSD) 
mechanism (see discussion in [|54l. l55|] ). The following weighted average value for this 
half-life is obtained: 

ri/2 = (7.1±0.4) 



2iLl2i[il]). In addition, here we have 

f2v\. First of all, the original result was 



In framework of High State Dominance (HSD) mechanism (see [|56l. l57|] ) the following 



average value can be obtained, T1/2 = (7.6 ±0.4) 



10^ 



IOOt 



Mo - i'^ORu (0+; 1130.29 keV) 

The transition to the 0+ excited state of ^^*^Ru was detected in five independent 
experiments. The results are ingood agreement, and the weighted average value for 
half-life is (using resuks from [[26,1271, |29m30(]): 



ri/2 = 5.9^0;? -lO-j. 



^20, 



The result from ||28|[ was not used here because we consider the result from [|29l] as final 
result of TUNL-ITEP experiment. 



116 



Cd 



There are three independent "positive" results that are in good agreement with each 
other when taking into account the corresponding error bars. Again, we use here the 
corrected result for the half-life value from Ref. [|33n . The original half-life value was 



^ We do not consider the result of Ref. ["52*1 because a possible high background contribution to the "effect" 
was not excluded in this experiment. 



decreased by ~ 25% (see remark in section 3.5). The weighted average value is (SSD 
mechanism): 

ri/2 = (2.8±0.2)-10iV 

If the HSD mechanism is realised for the case of ^ ^^Cd as well, then the adjusted half-life 
value is ri/2 = (3.0 ± 0.2) • lO^^ y. 



i^^Te and ^^^Te 

Long time there were only geochemical data for these isotopes. Although the half- 
life ratio for these isotopes has been obtained with good accuracy (~ 3%) [|35il . the 
absolute values for 7^/2 of the individual nuclei are different from one experiment to 



the next. One group of authors 1341 bSl . 159^1 gives 7^12 ~ 0.8 ■ 10^^ y for ^■^'^Te and 



ri/2 ~ 2 ■ 10 y for ^^'^Te, whereas another group [Il5l . l35n claims Ti /2 ~ (2.5 — 2.7) • 10 



^21 



y and ri/2 ~ 7.7 ■ 10 y, respectively. Furthermore, as a rule, experiments with "young" 

samples (~ 100 million years) result in half-life values of ^^^Te in the range of ~ 
(0.7 — 0.9) ■ 10^^ y, while for "old" samples (> 1 billion years), half-life values in 
the range of ~ (2.5 — 2.7) ■ 10^^ y have been produced. It was even assumed that the 
difference in half-life values could be connected with a variation of the weak interaction 
constant Gp with time [|50|] . 

One can estimate the absolute half-life values for ^^^Te and ^^^Te using only very 
well-known ratios from geochemical measurements and the "present" half-life value 
of ^^Se (see section 3.3). The first ratio is given by ri/2(^^^Te) /ri/2(^^^Te) = (3.52 ± 
0.11) ■ 10-4 [35J, while the second ratio is given by T^I2{^^Hq) /Tii2{^'^Sq) = 9.9 ±0.6. 
This latter value is the weighted average value from three experiments with minerals 



containing both elements (Te and Se): 7.3 ±0.9 [l60|], 12.5 ±0.9 lU and 10 ±2 116 10 . 
It is significant that the gas retention age problem has no effect on the half-life ratios. 
Now, using the "present" ^^Se half-life value ri/2 = (0.92 ±0.07) ■ 10^0 y and the value 

9.9 ±0.6 for the ri/2(^^*^Te)/ri/2(^^Se) ratio, one can obtain the half-life value for 
i30xe: 

ri/2 = (9±i)-io2V 

Using ri/2(i3^Te)/ri/2(i28jg) ^ (3.52±0.11) ■ 10"^ one can obtain the half-life 
value for ^^^Te: 

ri/2 = (2.5±0.3)-102V 
Recently it was argued that "short"half-lives are more likely to be correct [3^, [37]. Using 



different "young" mineral results half-life value was estimated as(9.0± 1.4) ■ 10^0 y [13611 
and(8.0± 1.1) • lO^o y Q for i^Otc and (2.41 ±0.39) ■ 10^^ y ^ and (2.3 ±0.3) • 10^4 
y [l37n for ^^^Te (corresponding to the observed ri/2(^^°Te)/ri/2(^^^Te) ratio). 



The first indication of a positive result for ^^*^Te in a direct experiment was obtained 
in [38]. More accurate and reliable value was obtaind recently in NEMO-3 experiment 
[I39.I . The results are in good agreement, and the waighted average value for half-life is 



ri/2 = (6.8±1.2)-1020y 



And now, using the 7^/2 (^'^'^Te) 77^/2 (^^^Te) ratio, one can obtane half-life value for 
i28Te, 

ri/2 = (2.0±0.3)-102V 

We recommend to use these last two results as most precise and reliable half-life values 
for ^■^'^Te and ^^^Te, respectively. 



150 



Nd 



The half-life value was measured in three independent experiments EHSil]. Using 
Equation 1, and three existing values one can obtain ri/2 = (8.2 ±0.5) ■ lO^^ y. Taking 
into account the fact that > ^ S = 1.89 (see Equation 2) we finally obtain: 



'1/2 



(8.2±0.9)-10iV 



150t 



Nd - i^^Sm (Oj^; 740.4 keV) 



+. 



There is only one positive result from a direct (counting) experiment 



'1/2 



[1.33^«i6(.^a^)^o:^3(^3'^0]■10^V 



0.27/ 



^20, 



Preliminary result of this work was published in [16211 



238 



u 



There is only one positive result from a radiochemical experiment [|43h 

ri/2 = (2.0±0.6)-102V 



130 



Ba (ECEC) 



There is only one positive result from a geochemical experiment ||44|1 : 

ri/2 = (2.2±0.5)-102V 

CONCLUSION 

In summary, all "positive" 2v/3/3 -decay results were analyzed and average values for 
half-lives were calculated. For the cases of ^^^Te and ^^'^Te, so-called "recommended" 
values have been proposed. We strongly recommend the use of these values as presently 
the most precise and reliable. In particular, the accurate experimental 2v/3/3 -decay rates 
can be used to adjust the most relevant parameter in the framework of the QRPA model, 
namely the strength of the particle-particle interaction (gpp). Once accomplished, these 
values can be used in NME calculations for neutrinoless double beta decay p63Ll64|l . 
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TABLE 1. Present, "positive" 2v^^ decay results. Here, N is the number of useful events, 
S/B is the signal-to-background ratio. *' For > 1.2 MeV. **' After correction (see text). 
***) For SSD mechanism. ****' In both peaks. 



Nucleus 


N 


Ti/2, y 


S/B 


Ref., year 


48Ca 


- 100 
5 

116 


[A3t\\{stat) ±\A{syst)] ■ IQ^'^ 

4.213-3.1019 
[4.4+04 (ifflf ) ±0.4(iysf) ' ^^^^ 

Average value: 4.4+JJ j IQi" 


1/5 
5/0 
6.8 


[4], 1996 
[5], 2000 
[6], 2008 


^"Ge 


-4000 

758 

132 

132 
-3000 
- 80000 


(0.9±0.1)-102i 
Ll^oj-lO^i 
0.931^-^ 10^1 
1.2+|]-2.102i 
(1.45 ±0.15)- 10^1 
[\.lA±Qm{stat)tl\l{syst)] ■ lO^^ 

Average value: (1.5 ±0.1) 10^^ 


-1/8 
-1/6 
^4 
^4 

- 1.5 

- 1.5 


[7], 1990 
[8], 1991 
[9J, 1991 
[10], 1994 
[11], 1999 
[12], 2003 


82Se 


149.1 
89.6 
2750 


[0.83 ±0.10(rfflf) ±0.07 (iyrf)] ' lO^** 

1.08lO-2«.1020 
[0.96 ± 0.03 (ifflf) ±0.1 (.yyif)] • 10^° 
(1.3 ±0.05) •102"(geochem.) 

Average value: (0.92 ± 0.07) lO^" 


2.3 

-8 
4 


[13], 1998 
[14], 1992 
[1], 2005 
[15], 1986 


96Zr 


26.7 
453 


[2.1+Jjj(ifaf)±0.2(iyrf)]-10''^ 

[2.35 ±0.14(sfa0 ±0.19(iyrf)] • 10^'* 
(3.9±0.9)-10''^(geochem.) 
(0.94 ± 0.32) •10''* (geochem.) 

Avprasjp valiip- + 2~1 . lO^' 


1.9*) 
1 


[16], 1999 
[a, 2009 
[17], 1993 
[18], 2001 


ioomo 


-500 

67 

1433 

175 

377 

800 
219000 


11 -6^0 8 -10'** 
[7.3 ± 0.35(rfflf ) ± 0.8(iyif)] • 10'*^**) 
7.6^^2.101^ 
[6.75;[J-^^(5fflf)±0.68(iyrf)] • 10'^ 
[7.2± l.l(5faf) ± 1.8(iyif)] • 10^** 
[7.11 ±0.02(,?fflf) ±0.54(iyif)] ' lO'^***' 
(2.1±0.3)-10'Vgeochem.) 

Average value: (7.1 ±0.4) 10^^ 


1/7 
7 
3 

1/2 
10 
1/9 
40 


[19], 1991 
[20], 1991 
[21], 1995 
[22], 1997 
[23], 1997 
[24], 2001 
[1], 2005 
[25], 2004 



TABLE 1. continued. 



ioomo - 
iooRu (0+) 


133****^ 
153****^ 
19.5 
35.5 
37.5 


6.i;|:?-io2o 

[9. 3tl ^{st at) ± 1. 4{sy St)] -10^° 
[5.9t[:](5fflf)±0.6(iyrf)]-102'' 
\5.5+\iistat)±0.3(syst)] ■ lO^" 
[5.7tol{stat) ±0.8isyst)] ■ lO^" 

Average value: 5.9+JJ^ 10^" 


1/7 
1/4 
~8 
~8 
-3 


[26], 1995 
[27], 1999 
[28], 2009 
[29], 2009 
[30], 2007 




- 180 
9850 
174.6 
1370 


2.6l«:^10i9 
[2.9±0.06{stat)toi{syst)] ■ lO'*^ 
[2.9±0.3{stat)±0.2{syst)] ■ lO'*^**) 
[2.8 ±0.1 (ifflf) ±0.3(iyif)] ■ 10''^***) 

Average value: (2.8 ± 0.2) IQi" 


^1/4 
-3 
3 

7.5 


[31], 1995 
[32], 2003 
[33], 1996 
[6], 2008 






- 2.2 • 10^4 (geochem.) 
(7.7±0.4)-102'* (geochem.) 
(2.41 ±0.39) ■ 10^4 (geochem.) 
(2.3 ±0.3)- 10^4 (geochem.) 

Recommended value: (2.0 ± 0.3) 10^'* 




[34], 1991 

[35], 1993 
[36], 2008 
[37], 2008 


i30xe 


260 
236 


[6.1 ± lA{stat)tl-lisyst)] ■ 10^^ 

[6.9±0.9{stat)tljisyst)] ■ 10^° 
~ 8 • lO^" (geochem.) 
(27 ±1)- 10^0 (geochem.) 
(9.0±1.4)-102'' (geochem.) 
(8.0 ±1.1)- 10^0 (geochem.) 

Recommended value: (6.8 ±1.2) lO^" 


1/8 
1/3 


[38], 2003 

[39], 2009 
[34], 1991 
[35], 1993 
[36], 2008 
[37], 2008 


150Nd 


23 
414 
2018 


[I8.8tti{stat) ±l.9{syst)] ■ 10^^ 
[6J5toAl{stat) ±0.68{syst)] ■ lO'** 
[9.11 toil {stat ) ± 0.63 (syst)] ■ 1 ' ^ 

Average value: (8.2 ± 0.9) lO^^ 


1.8 

6 

2.8 


[40], 1995 
[23], 1997 
[41], 2008 


150Nd - 

i^^Sm (0+) 




[l33t'of,{stat)tlT3{^yst)]-\Q^^ 
Average value: 1 33+JJ lO^" 


1/5 


[42], 2009 


238u 




(2.0 ± 0.6) • 10^1 (radiochem.) 




[43], 1991 


ECEC(2v) 




(2.2 ±0.5) -10^1 (geochem.) 




[44], 2001 



